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Abstract. We report wide class of exact solutions of the modified Gross-Pitaevskii equation 
(GPE) in 'smart' Jacobi elliptic potentials: V(^) = -Vosn(^,m), V(^) = -Vocn(^,m), and 
y(f) = -VodnCf, m) in the presence of external source. Solitonlike solutions, singular 
solutions, and periodic solutions are found using a recently developed fractional transform: 

p(f ) = ^^^^yr , where / is the respective Jacobi elliptic function and the amplitude parameters 
A, B, and D nonzero. These results generalize those contained in (Paul T, Richter K and 
Schlagheck P 2005 Phys. Rev. Lett. 94, 020404) for nonzero trapping potential. 
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1. Introduction 



In a mean-field approximation, the dynamics of a dilute-gas Bose-Einstein condensate (BEC) 
can be captured by the cubic nonlinear Schrodinger equation (NLSE) with a trapping potential 
|[Tl|2l|3l. The various traps which are used to contain the BEC have spurred the solutions of 
the NLSE with new potentials |i4j We consider the mean-field model of a quasi-one- 
dimensional BEC trapped in a 'smart' potential in the presence of an external source||6J 
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1^ + KexpiixiO - i0t), (1) 



where if/{x, t) represents the macroscopic wave function of the condensate and V{^) is an 
experimentally generated macroscopic potential. The parameter g indicates the strength of 
atom-atom interactions and it alone decides whether Eq. (1) is attractive (g = -1, focussing 
nonlinearity) or repulsive (g - 1, defocussing nonlinearity). Here, K and are real constants 
related to the source amplitude and the chemical potential, respectively, and x(0 is a real 
function of ^ = a(x - vt), a and v being two real parameters. In the field of nonlinear optics, 
Eq. (1) may describe the evolution of the local amplitude of an electromagnetic wave in the 
spatial domain, in a two-dimensional waveguide where t becomes the propagation distance 
and X becomes the retarded time, and the system is driven by an external plane pump wave. 
The Jacobi elliptic potential may describe a transverse modulation of the refractive index in 
the waveguide. 

As is well-known, Eq. (1) is not integrable if KVq + 0, and only small classes of explicit 
solutions can most likely exist. For y(^) - 0, Eq. (1) is a cubic NLSE with a source, and 
exact rational solutions using a fractional transform are found in Ref. |7|. And in Ref. ifTTIl 
periodic solutions without source have been reported. More recently, in Ref.|81, a class of 
exact solutions of Eq. (1) for y(^) - -Vosn^C^, m) have been reported. In particular, the 
rational solutions of the fractional transform: p(^) = f^i/^' where / = sn have been reported 
for B = 0. This is due to the form of the trapping potential. Nonetheless, in the present paper 
we find rational solutions of the type p(^) - where / being the respective Jacobi elliptic 

function with all the amplitude parameters A, B, and D nonzero. These results generalize 
those contained in Ref |6| for nonzero trapping potential. The choice of a smart potential 
y(^) allows one to construct a large class of exact solutions, as done in a number of works for 
the cubic GP equations jO] \lU\ [TT| . In the present work, we consider three different potentials 
in the GP equation: V(^) = -yosn(^,m), V{^) = -yocn(^,m), and V{^) = -yodn(^,m) 
in the presence of external source, and find exact travelling wave solutions of Eq. (1) with 
K 0. The choice of these three different 'smart' potentials is motivated by the following 
facts. Firstly, the potential V{^) = -yosn(^, m) in the limit m ^ is V(^) = -yosin(^) 
which is similar to the standard optical lattice potential |12, 13 1. Secondly, the choice of 
the potential V(^) = -y()dn(^, m) in GP equation mimics| 14| the harmonic potential that was 
used to achieve BEC experimentally. The third potential, we hope it is relevant to the available 
experimental conditions to achieve BEC. 



2. EXACT SOLUTIONS OF THE GPE IN 'SMART' PERIODIC POTENTIAL WITH 
SOURCE 

The travelling wave solutions of Eq. (1) with potential V{^) are taken to be of the form 
i//(x, t) - p(^)e''*^['^(^""')l""*'. Inserting this expression for 4r(x, t) in Eq. (1) and separating the 
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real and imaginary parts, and integrating the imaginary part, one gets 

V C 

+ — (2) 

where C is a constant of integration. In order that the external phase be independent of 41, we 
consider only solutions with C = 0, to obtain 

P + ( :j—)P --TP 2~P ~~2 

Below we consider three different 'smart' Jacobi elliptic potentials |151 and find exact 
solutions. 

Case(l):- V(^) = -yosn(^,m). Eq. (3) reads as 



,P' 

a 

Substituting 



y + 2g 3 2Vo 2K ^ 

( —)P - 4p + ^snp - ^ = 0. (4) 



p(^)=Ai +Bisn(^,m) (5) 

in Eq.(4) and equating the coefficients of equal powers of sn(^, w) result in relations among 
the solution parameters A 1 , and B\, and the equation parameters Vo, g, K, a, and 0. We find 
that 

a = -;— — , (6) 
1 + m 



V() I ma^ 

Ai = T-^, Bi = J . (7) 

3a^lgm ^ g 

From Eq. (7) it follows that Vq > and gf > implying the GPE with repulsive nonlinearity. 
Case(II):- V{^) = -yocn(^,m). Eq. (3) reads as 

„ ,v^ + 20^ 2g 3 2Vq 2K ^ 
P + ( - -|p^ + ^cnp - ^ = 0. (8) 

Substituting 

p(^)=A2 + B2cn(^,m) (9) 

in Eq.(8) and equating the coefficients of equal powers of scn(^, m) result in relations among 
the solution parameters A 12, and Bi, and the equation parameters Vq, g, K, a, and 0. We find 
that 

20 + 

" - 1^' '''^ 



5:=J^. (11) 



3a -\J-gm V g 

From the positivity of we conclude that cn solutions exist for Vq > and g < Q. The 
condition gf < corresponds to the GPE with attractive nonlinearity. 
Case(lll):- V(^) = -yodn(^,OT). Eq. (3) reads as 

„ ,v^ + 20^ 2g 3 2Vo^ 2K ^ 
p" + ( ^)p - ^p' + ^dnp - — = 0. (12) 

Substituting 

p(^) = A3 + B3dn(^, m) (13) 
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in Eq.(12) and equating the coefficients of equal powers of dn(^, m) result in relations among 
the solution parameters A3, and B3, and the equation parameters Vo. 9^ K, a, and 0. We find 
that 

2 20 + t;2 

or = -, (14) 

m-2 



A3 = -^, ^3=J^. (15) 
Here we conclude that dn solutions exist only for Vq > and g <0. 

3. Rational solutions 

In order to obtain Lorentzian-type of solutions of Eq. (3) we use a fractional transform 

where / is the respective Jacobi elliptic functions. Again we obtain the Lorentzian-type of 
solutions of Eq. (3) for three different 'smart' potentials. 
Case(I):-y(^) - -yosn(^,m). Eq. (3) reads as 

„ y + 20 2g 3 2^0 2K 
P + ( T—)P - -jP + —r^^P 2" = 0- (^^) 

Substituting 

= ^^iP^ (18) 
1 + Z)isn(^,m) 

in Eq.(17) and equating the coefficients of equal powers of sn(^,m) will yield the following 
consistency conditions. 

2B4 + 2AaD\ + YAa - - ^ = 0, (19) 

a'- 

2mBA.D\ - = 0, (20) 



a'- 

dmBiPx + ^B4£>i = 0, (21) 



6mB4 + B40i(r - m - 1) - ^A^bI + ^B^D^ = 0, (22) 
-2mAiDi + B^DxiTC - 3w - 3) + ^A^D] + ^B - = 0, (23) 

-454(1 +m)+ A^DjiF - m - 1) + m - -^aIb^ + -^A^Di - —D\ ^) 

Or^ Or^ Or^ 

A^Di{2T + m + 1) + — 5^A4 - -^-Di = 0. (25) 

a a 

From the above consistency conditions we obtain the following relations. 

^4 = ^, (26) 

3m{m + \)a^ + dmTa^ — ^ 



Substituting 
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where T = Here, we would like to emphasize that these results generalize those 

contained in Ref IQ, for nonzero trapping potential. This stems from the fact that the constant 
B in expression (18) is nonzero, which follows from the choice of our 'smart' potential in Eq. 
(1). 

Case(II):-y(f) = -yocn(^,m). Eq. (3) reads as 

p" + ( r^)p - -|P + —T'^'^P - — = 0. (28) 

Pi^)-- T.T}r. (29) 
1 + D2cn(^, m) 

in Eq.(28) and equating the coefficients of equal powers of cn(^, m) will yield the following 
consistency conditions. 

2^5(1 -m) + 2A5D?(1 -m) + YAs - ^AI - ^ = 0, (30) 

-2mB45Dl-^Bl^Q, (31) 

-bmBiD2 + ^B^dI = 0, (32) 

-6wB5 + BsDi,{Y + 2m - 1) - ^A^bI + ^B5£>2 = 0, (33) 

2mAiD2 + BiD2(2T + 6m - 3) + —^A^dI + -^Bs - —rOl = 0, (34) 

, 6fl , 4y() , 

-4B5(1 - 2m) + AiD\{Y + 2m - 1) + TB, - -^A\B^ + — f A5D2 - —rD\ #39) 

ff- a'^ 

A,D2(2r - 2m + 1) + - = 0. (36) 

From the above consistency conditions we obtain the following relations 

\%mK 



2V 

3m( l - 2m)a^ + 6mra^ H — t- 



(37) 



3a^m j m 3ma^ 
B5 = J , £>2 = ■ (38) 



Case(III):- y(^) = -yodn(f , m). Eq. (3) reads as 



,".,fi±a»),_?£p>.?}^a„p-?£ = o. (39) 



Substituting 



1 + D3dn(^, m) 

in Eq.(39) and equating the coefficients of equal powers of dn(^, m) will yield the following 
consistency conditions. 

2A6D2(m - 1) + rA6 - ^A^ - ^ = 0, (41) 

-2BeDl -^Bl^O, (42) 
2Vo , 

-4B6D3 + —^B(,dI = 0, (43) 
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Figure 1. Singular solitary wave solution for a = \ and Vo = I, and 3=1. 



-4B6 - 2B6£>3 + BeDliY - m + 2) - %AeBl + ^B^D^ = 0, (44) 

lAeDi + B^D^ilT - 3m + 4) + ^AgD^ + - = 0, (45) 

a'^ a 

-2Bf,(m - 2) + IBeDi + AeOjiT -m + 2)+ TBg - ^A^Bg + ^AgDs - ^Dj 4«) 

A6D3(2r + m - 2) + - = 0. (47) 

From the above consistency conditions we obtain the following relations 

6K 

Aa = ^, (48) 

Q'2(2r + m - 2) + ^ 



2ck2 q.2 2a^ 



3.1. Trigonometric solutions 

From the consistency conditions that arise from the first two 'smart' periodic potentials, we 
conclude that the limit m = is forbidden as the amplitude parameters A, B, and D will be 
zero. On the other hand, for V{^) - -yodn(f , m) case, only flat background solutions wiU be 
possible for m - limit. 

3.2. Solitonlike solutions 

Here, in this subsection, we describe the solitonlike solutions that are obtained from the 
solutions in sn(^, m) and cn(^, m) in the limit m - 1, in detail. In the limit m = 1, V{^) 
becomes an array of well separated kink-type of potential barriers: V{^) = -yotanh(^). Then 
we have the following relations 

18/: 
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Figure 2. Non-singular solitary wave solution for a = l,yo = -6, ^ = -1 and K = 1/2. 



And the strength of the source is 
K = 



yo[6(r + - 2y2] (63(T - 2)a* - 2y2^ 



1O8q'30 



1/2 







y2 



As a special case, if we set ff = 1 and yy = 1, then we get A4 — ^ffj' ^4 - 3 / -\/g, and 
Di - -3. This results in a solitonlike solution 

'^^^'^ l-3tanh(^) ^ ' 

This set corresponds to the singular solution for repulsive case i.e., g > 0. The singularity 
of the pulse profile may correspond to the beam power exceeding material breakdown due 
to self-focussing lfT6l [TtI [TSl [191 . Figure (1) depicts a surface plot of this solution for the 
parameter values given in the figure caption. 

Another interesting solitonlike solution is obtained from the solution in cn(^, w) for 
m - 1. In the limit m - I, V(^) becomes an array of well separated secant hyperbolic 
potential barriers: V(^) = -yosech(^).Then we have the following relations 

A, = (53) 

3(2r-l)a2 + 9i 

Bs^'—^T^g, D2^^. (54) 

As a special case, if we set a = 1, Vo - -6, g - -I and K = 1/2 then we get A5 = 
£5 - -(1/2), and Di - -(1/2). This results in a solitonlike solution 

As + fi5sech^(^) 
1 - 0.5sech(^) 

This set corresponds to the non-singular solution for attractive case i.e., g <Q. The same has 
been depicted in Fig. 2 for the parameter values given in in the figure caption. 



6r+69 ' 



4. Conclusions 



In conclusion, we have shown the existence of wide class of exact solutions for the modified 
GP equation in 'smart' periodic potentials with an external source. The Lorentzian-type of 
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solitons are obtained with the aid of a fractional transform. Our analysis applies to both 
attractive and repulsive cases of GP equation. Furthermore, our rational solutions generalize 
those contained in Ref |6|, for nonzero trapping potential, because of our choice of 'smart' 
potential. We hope that these potentials may be experimentally realizable, to achieve BEC. 
Although not presented here, the stability of these wide class of solutions can be checked using 
a semi-implicit Crank-Nicholson finite difference method |7|, as the much used numerical 
techniques based on fast Fourier transform (FFT) requires the FFT of the source, which is 
costly. 
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